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Publications 

This  contract  was  awarded  for  research  into  the  development  of<*:,'^' 
numerical  and  analytic  techniques  useful  in  solving  problems  arising 
in  continuum  mechanics,  especially  fluid  mechanics,  and  particularly 
in  areas  in  which  the  PI  and  co-PIs  have  special  expertise.  Projects 
that  were  pursued  during  tenure  of  the  award  including  the 
following. 


I.  Dynamical  Processes  in  Freezing  Materials.  Two  main  areas  were 
attacked  during  the  period  of  the  contract:  (a)  the  development  of 
continuum  theories  of  freezing  processes,  and  (b)  the  solution  of 
moving  boundary  problems  connected  with  the  freezing  of  materials, 
(a)  Continuum  Theories.  A  system  near  the  meliing/freezing  point 
typically  contains  matter  in  three  forms;  solid,  liquid  and  mush,  the 
last  being  a  mixture  of  solid  and  liquid  phases.  As  a  system  evolves, 
it  converts  each  of  these  forms  into  one  of  the  others,  i.e.  the  solid- 
liquid,  solid-mush  and  liquid-mush  fronts  all  move  relative  to  the 
material.  The  objective  is  to  develop  a  physically  well-founded 
theory  that  will  correctly  describe  the  evolution  of  each  domain  and 
the  motion  of  the  interfaces  separating  them  (see  Hills  &  Roberts, 
1987a,b  1988a,b,c;  Loper  &  Roberts,  1987a,b;  Roberts  &  Loper  1988). 
Complications  that  arise  include  coarsening  (the  increase  in  the  mean 
size  of  the  solid  grains  as  the  mush  ages;  see  Hills  &  Roberts,  1990), 
coring  (the  fact  that  the  material  currently  freezing  onto  the  grains  of 
an  alloy  does  not  have  the  same  composition  as  that  deeper  in  the 
grain,  which  froze  earlier;  see  Hills  &  Roberts,  1988b,c),  effects  of 
surface  tension  (particularly  the  freezing  temperature,  which  is 
depressed),  and  many  other  processes.  Theories  were  developed 
that  included  the  effects  just  specifically  mentioned.  These  were 
tested  by  analytic  solutions  to  simple  situations,  and  the  results  gave 
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confidence  that  the  theories  were  physically  robust.  Related  ideas 

were  applied  to  other  mixture  theories  (Hills  &  Roberts,  1987c,  1989) 
(b)  Specific  Solutions.  To  determine  the  evolution  of  a  system  near 
the  freezing  point,  requires  that  complicated  Stefan  problems  be 

solved  to  determine  the  location  and  motion  of  the  interfaces 
mentioned  in  (a)  above.  One  challenging  problem,  to  which 

previously  only  solutions  involving  such  physical  absurdities  as 
superheated  ice  had  been  proposed,  was  solved  at  UCLA  during  the 
period  of  the  contract  (Hills  &  Roberts,  1991);  it  concerns  the 

evolution  of  system  cooled  below  the  freezing  point  on  its  boundary 
but  healed  internally  by  (e.g.  radioactive)  sources. 

2.  Nonlinear  Waves.  A  number  of  projects  were  undertaken 
concerning  the  interaction  and  evolution  of  a  sea  of  small  amplitude 
waves,  A  Fokker-Planck  equation  was  developed  that  generalizes  all 
known  previous  Focker-Planck  equations,  and  with  the  help  of  which 
the  frequency  renormalization  created  by  nonlinear  interactions 
could  be  explicitly  calculated  in  a  first  approximation  (Putterman  & 
Roberts,  1988).  The  particular  case  of  acoustic  turbulence  was 
studied  both  for  classical  systems  (so  generalizing  well-known 
Russian  work)  and  quantum  systems  (so  shedding  more  light  onto 
the  interface  between  classical  and  quantum  physics).  Large-scale 
numerical  simulations  were  also  undertaken  (Haeri  et  al.  1992). 
Further  work  on  ocean  waves  revealed  the  existence  of  a  second 
sound  phenomenon  (Larraza  1987,  Larraza  et  al.  1990). 

Another  project  that  still  engages  the  attention  of  the  PI 
includes  three-dimensional  integrations  of  the  non-linear 
Schrodinger  equation  as  a  simple  model  of  superfluid  helium  (the 
Bose  condensate  model)  ;  see  Jones  et  al.  (1986),  Jones  &  Roberts 
(1987). 


3.  Convection.  A  number  of  investigations  of  MHD  stability  were 
completed.  Magnetoconvection  was  studied,  both  for  compressible 
and  Boussinesq  fluids  (Jones  &  Roberts,  1990;  Roberts,  1991).  The 
well-known  studies  of  the  resistive  instabilities  of  plasma  equilibria  . 
was  generalized  to  highly  rotating  systems  (Kuang  &  Roberts  1990,  . 
1991,  1992).  A  nonlinear  study  of  the  arotropic  annulus  (a  crude  i 
model  of  the  atmosphere  of  Jupiter)  was  undertaken:  see  Lin  (1990). 
Mulder  (1989)  developed  an  new  multigrid  approach  for  the  solution 
of  convection  problems. 
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collaboration  with  Robert  Miura 

Uuiv,,  Los  An^'cies.  . 

Per  Telecon  ONR  6f2hl93  cl 


Lange  continued  his  long  term 
on  perturbation  methods,  and 


1^ 


Avail  aiuljor 
i>pi..c  f  a! 


<4 


published  an  investigation  of  the  asymptotics  of  differential- 
difference  equations  (Lange  &  Miura,  1991).  Lange  has  also 
continued  his  collaboration  with  Don  Smith  on  integral  equations 
involving  a  small  parameter.  Together,  they  have  completed  a 
lengthy  study  of  the  application  of  singular  perturbation  methods  to 
such  problems  (Lange  &  Smith,  1988).  Unusual  kind  difficulties  arise 
from  the  theory  of  elasticity  of  elastic  membranes  and  shells.  These 
have  significant  technological  applications  to  tubular  chemical 
reactors.  Lange  &  Weinitschke  (1991a,  b)  have  developed  the 
necessary  mathematical  apparatus  to  solve  these  problems. 

5.  Vorticity.  During  tenure  of  the  award,  Anderson  continued  his 
investigations  into  the  creation,  interaction  and  merging  of  vorticity. 
The  results  of  such  study  are  technilogically  significant,  since 
vorticity  creation  and  shedding  has  a  significant  bearing  on  the  lift 
and  drag  of  a  moving  body.  It  can  also  be  argued  that  the  vorticity 
contains  the  physics,  and  that  the  primary  task  in  integrating  the 
equations  of  motion  is  that  of  solving  the  vorticity  equation,  and 
subsequently  deriving  the  flow  subsequently  from  the  relevant 
Poisson  equation,  perhaps  by  domain  decomposition  methods 
(Anderson  1988a).  The  direct  solution  of  the  vorticity  encounters 
difficulties  in  satisfying  the  boundary  conditions.  Anderson  has 
developed  new  techniques  for  overcoming  these  (Anderson  1988b, 
1989a),  and  has  applied  the  resulting  methods  to  problems  involving 
vortex  shedding  (Anderson  et  al.  1990;  see  also  Anderson  et  ai, 
1988a)  as  well  as  to  the  interaction  of  vortex  rings  (Anderson  et  al., 
1988b;  see  also  Anderson,  1989b).  Anderson  (1990)  also  made  a 
number  of  suggestions  concerning  fast  particle  simulations,  as  well  as 
alternative  implementations  (Anderson,  1992)  of  the  multipole 
method  of  Greengard  &  Rokhlin  (1991).  Cottet  (1990a,  b)  also  added 
to  the  theory  of  deterministic  vortex  methods  for  the  Navier-Stokes 
equations,  with  particular  reference  to  implementation  of  the 
boundary  conditions,  ’hrough  a  particle-grid  superposition  method. 

6.  Shocks.  During  tenure  of  the  award,  many  applications  have  been 
made  by  Osher  and  his  collaborators  of  essentially  non-oscillatory 
schemes  for  shock  capture:  Fatemi  et  al.  (1991a,b),  Harten  (1988, 
1990),  Harten  &  Chakravarthy  (1992),  Harten  et  al.  (1987),  Mao 
(1989a,b,  1990),  Marquina  (1992),  Osher  &  Shu  (1989,  1991)',  Shu  & 
Osher  (1989),  Shu  et  al.  (1992),  Yang  (1990).  These  can  be  used  to 
compute  the  motion  of  shocks  whose  speed  depends  on  their 
curvature,  for  which  a  Hamilton-Jacobi  formulation  is  particularly 
convenient  (Bardi  &  Osher,  1991;  Osher,  1989,  1991;  Osher  &  Shu, 


1991).  This  has  applications  to  solving  problems  involving  crystal 
growth  and  those  involving  the  propagation  of  flames  or  detonations 
(Ton  et  al.,  1992).  The  general  problem  of  calculating  the  motion  of 
an  interface  is  a  challenging  one;  other  contributions  have  been  made 
by  Mao  (1991a,b),  Mulder  ei  al.  (1991).  Chu  &  Wu  (1991)  developed 
the  theory  of  magnetohydrodynamic  shocks  further.  Engquist  et  al. 
(1989b)  developed  nonlinear  numerical  filters  for  shock  calculations. 
Mostrel  (1988)  developed  new  numerical  schemes  for  transonic  flow 
problems.  Sjogreen  (1988)  reported  on  shock  calculations  using 
algorithms  developed  for  parallel  computation.  Engquist  &  Huynh 
(1989)  developed  iterative  pre-conditioned  Newton-type  methods 
for  shock  calculations. 

Several  investigations  were  made  of  problems  encountered 
when  conservation  laws  are  solved  computationally.  Lafon  &  Osher 
(1991)  explored  high  order  filtering  methods  for  the  solution  of  high 
order  systems  of  hyperbolic  conservation  laws.  Donat  &  Osher 
(1990)  investigated  the  propagation  of  errors  in  nonlinear 
approximations  to  conservation  laws.  Marquina  (1990)  used  a  local 
piecewise  hyperbolic  reconstruction  of  numerical  fluxes  in  the 
solution  of  nonlinear  conservation  equations.  Durlofsky  et  al.  (1990) 
developed  a  triangle-based  TVD  scheme  for  the  solution  of 
hyperbolic  conservation  laws.  The  incorporation  of  a  maximum 
principle  into  a  triangular  stencil  was  implemented  by  Liu  (1992). 

Engquist  &  Sjogreen  (1990)  gave  new  methods  for  solving  stiff 
hyperbolic  systems.  Brenier  (1991)  added  to  the  theory  of 
hyperbolic  differential  equations  that  have  oscillatory  coefficients. 
Papageorgiou  &  Smyrlis  (1991)  studied  the  route  to  chaos  for  the 
Kuramoto-Sivashinsky  equation. 

7.  Turbulence .  Henshaw  has  almost  completed  a  fourth  order 
accurate  difference  code  for  Navier-Stokes  equations,  using 
overlapping  grid  techniques.  One  can  use  rather  general  boundary 
conditions.  Adaptivity  will  be  added.  As  part  of  her  thesis  Jennifer 
Hou  (nee  Chang)  has  compared  finite  differences  and  spectral 
methods  for  model  problems.  She  has  also  implemented  the  methods 
on  the  connection  machine.  Henshaw  and  Kreiss  have  made 
perturbation  calculations  for  2D  turbulent  flows.  The  following 
results  arv  interesting; 

(a)  The  effect  of  the  perturbations  depend  very  much  on  the 

frequency  of  the  perturbation; 

(b)  Fully  turbulent  flows  are  stable  against  perturbations  in 

short  time  intervals; 


(c)  In  long  time  intervals  the  effect  of  a  change  in  Reynolds 
number  is  large.  For  example,  changing  R  fromlO^  to  2.10^ 
changes  the  large  scale  by  0(1).  To  be  able  to  make  the 
calculations  more  efficient  D.  Brown  has  implemented  our 
code  on  the  connection  machine  at  Los  Alamos  National 
Labratory. 

8.  Systems  having  Two  Widely  Different  Time  Scales.  Kreiss  & 
Lorenz  (1992)  have  completed  a  paper  on  slow  solutions  for  highly 
oscillatory  problems".  They  succeeded  in  connecting  the  bounded 
derivative  principle  with  asymptotic  expansions  and  the  existence  of 
a  slow  manifolds.  Also,  sharp  estimates  of  the  effect  of  the  fast  scale 
on  the  slow  time  scale  were  obtained.  Kreiss  (1992)  has  completed  a 
large  survey  article  on  problems  with  different  time  scales.  It 
discusses  the  connection  between  ordinary  and  partial  differential 
equations  and  also  treats  applications  in  Geophysics  and  Plasma 
physics.  Kreiss  et  al.  (1991)  have  presented  a  new  discussion  of 
convergence  of  the  solutions  of  the  equations  governing  a 
compressible  fluid  to  the  solutions  of  the  incompressible  Navier- 
Stokes  equations.  Browning  et  al.  (1992)  have  made  a  comparison  of 
the  differential  systems  and  the  numerical  methods  used  in  the 
computation  of  smooth  oceanographic  flows.  For  highly  oscillatory 
problems  the  time  integrator  is  still  an  obstacle.  Some  of  these 
difficulties  are  discussed  in  a  forthcoming  paper  by  Kreiss  & 
Abrah.''.msson  (1992).  A  paper  by  Kreiss  &  Scherer  (1992)  on  the 
method  of  lines  for  hyperbolic  differential  equations  wil  appear 
shortly.  B.  Gustafsson,  J.  Oliger  &  H.-O.  Kreiss  have  started  work  on  a 
book  on  Numerical  Methods  for  Time  Dependent  Problems. 

9.  Initial  boundary  value  problems.  Lixin  Wu  has  defended  his 
thesis:  "Stability  of  difference  approximations  for  initial  boundary 
value  problems."  Kreiss  has  worked  on  a  new  definition  of 
wellposedness  for  initial  boundary  value  problems.  It  is  based  on  the 
resolvent  condition  and  it  is  the  weakest  definition,  which  is  stable 
against  lower  order  perturbations.  Also,  the  discrete  analog  of  this 
definition  is  very  effective  for  the  method  of  lines. 

10.  Far-fie Id  Boundary  Conditions.  When  the  solution  of  hyperbolic 
equations  are  required  in  infinite  domains,  it  is  necessary  in  practice 
to  cut  off  the  domain  at  great  distance  and  to  add  boundary 
conditionf  at  the  cut-off  that  preserve  the  essential  condition  that 
the  waves  are  outgoing.  Reflection  in  the  far  field  would  lead  to  a 
return  of  the  waves  to  the  source  regions  and  to  the  corruption  of  the 


solution.  The  implementation  of  far-field  boundary  conditions  had 
been  of  special  interest  to  Engquist  and  his  collaborators  for  several 
years,  and  further  progress  was  made  during  the  period  of  this 
award.  The  main  contributions  were  applications  to  fluid  dynamics 
and  steady  problems.  See  Engquist  (1989),  Engquist  &  Gustafsson 
(1987)  and  Engquist  &  Halpern  (1988,  1990).  Global  boundary 
conditions  for  use  with  fast  Helmholtz  solvers  were  studied  by 
Engquist  et  al.  (1989a). 

11.  Other  Computationally-Oriented  Work.  Bamberger  et  al. 
(1988a,b)  developed  methods  for  solving  parabolic  equations  in 
heterogeneous  media,  Rizzi  &  Engquist  (1987)  surveyed  the  theory 
and  practice  of  computaional  fluid  mechanics,  concentrating  mainly 
on  specially  selected  topics. 

Conferences 

The  UCLA  Workshop  on  Vortex  Methods  was  held  at  UCLA 
during  May  20  through  May  22,  1987  and  was  supported  by  this 
Contract.  It  was  organized  by  Chris  Anderson  and  brought  together 
people  investigating,  by  computational  and  analytic  methods,  the 
evolution  of  vortical  structures.  It  resulted  in  a  book  in  the  Springer 
series  Lecture  Notes  in  Mathematics;  see  Anderson  &  Greengard 
(1988). 

The  Workshop  on  the  Nonlinear  Dynamics  of  Rotating  Magnetic 
Systems  was  held  at  UCLA  during  August  4  through  8,  1987  and  was 
supported  by  this  Contract.  It  was  organized  by  Paul  Roberts  and 
brought  together  41  people  interested  in  this  interesting  branch  of 
MHD,  one  that  has  wide  significance  in  planetary  and  extra¬ 
terrestrial  magnetism.  The  report  on  the  meeting  (Roberts,  1987) 
was  widely  distributed. 

The  Second  International  Symposium  on  Domain  Decomposition 
Methods  was  held  in  UCLA  during  January  14  through  16,  1988  and 
was  partially  supported  by  this  Contract.  Tony  Chan  was  the 
organizer  on  site  and  one  of  the  editors  of  the  final  published 
proceedings  (Chan  et  al,  1989).  Nearly  100  people  attended  the 
meeting,  and  40  featured  lectures  presented  state-of-the-art 
reviews,  covering  theory,  algorithms,  parallel  implementations,  and 
applications. 


References^ 

Anderson,  C.R.  1988a  In  Domain  Decomposition  Methods,  (Eds  T. 
Chan,  R.  Glowinski,  J.  Periaux  &  O.  Widlund).  pp]29-139.  SIAM 
Publications. 

Anderson  C.R.  1988b  In  Mathematical  Aspects  of  Vortex  Dynamics. 
(Ed.  R.F.  Caflisch).  ppl44-160.  SIAM  Publications. 

Anderson,  C.R.  1989a  J.  Comp.  Phys.  80,  12-91 . 

Anderson,  C.R.  1989b  Comm.  Pure  Appl.  Math.,  42,  1123-1139. 
Anderson,  C.R.  1990  In  Computational  Methods  in  Applied  Science 
&.  Engineering,  (Eds.  R.  Glowinski  &  N.  Lichnewsky).  ppl23-135. 
SIAM  Publications. 

Anderson,  C.R.  1992  "An  implementation  of  the  fast  multipole 
method  without  multipoles",  SISSC,  to  appear. 

Anderson,  C.R.  &  Greengard,  C.  (Editors)  1988a  Vortex  Methods, 
Springer  Lecture  Notes  in  Mathematics,  1387.  Heidelberg;  Springer. 
Anderson,  C.R.  &  Greengard,  C.  1988b  In  Vortex  Methods,  (Ed. 
Anderson,  C.R.  &  C.  Greengard).  pp23-37.  Springer. 

Anderson,  C.R.  Greengard,  C.,  Greengard,  L.  &  Rokhlin,  V.  1990 
Phys.  Fluids  A2,  883-885. 

Anderson,  C.R.  Greengard,  C.  &  Henderson,  M.  1988b  In  Vortex 
Methods,  (Ed.  C.R.  Anderson  &  C.  Greengard).  pp42-55,  Springer. 
Bamberger,  A.,  Engquist,  B.,  Halpern,  L.  &  Joly,  P.  1988a  SIAM  J . 
Appl.  Math.,  48,  99-128. 

Bamberger,  A.,  Engquist,  B.,  Halpern,  L.  &  Joly,  P.  1988b  SIAM  J. 
Appl.  Math.,  48,  129-154. 

Bardi,  M.  &  Osher,  S.J.  1991,  SIAM  J.  Math.  Anal.  22,  344-351. 
Brenan,  K.  &  Engquist,  B.  1988  Math.  Comp.  51,  659-676. 

Brenier,  Y.  1991  In  Proceedings  of  the  Uppsala  Conference  on 
Hyperbolic  Systems  of  Conservation  Laws,  (Ed.  B.  Engquist). 
Browning,  G.  Holland,  W.  &  Kreiss,  H.-O.  1992  "A  comparison  of 
differential  systems  and  numerical  methods  for  the  computation  of 
smooth  oceanographic  flows",  to  appear. 

Chan,  T.F.,  Glowinski,  R.,  Periaux,  J.  &  Widlund,  0.  (Editors)  1988 
Domain  Decomposition  Methods.  SIAM  Publications. 

Chu,  C.K.  &  Wu,  C.C.  1990  "Magnetohydrodynamic  shock  waves 
revisited".  Unpublished  CAM  Report  90-27 


t  Authors  supported  by  the  Grant  (which  include  the  Principal  Investigator, 
the  Co-principal  investigators,  other  faculty  members,  visitors  and  students) 
are  shown  in  bold  face  type.  The  phrase  "CAM  Report"  refers  to  the  Scries  of 
reports  on  Computational  and  Applied  Mathematics  produced  by  the 
Department  of  Mathematics  at  UCLA. 


Cottet  G.-H.  1989a  "Boundary  conditions  and  deterministic  vortex 
methods  for  the  Navier-Stokes  equations".  Unpublished  CAM  Report 
88-42. 

Cottet  G.'H.  1989b  "A  particle-grid  superposition  method  for  the 
Navier-Stokes  equations",  Unpublished  CAM  Report  89-19. 

Donat,  R.  &  Osher,  S.J.  1990  Comp.  Meth.  Appl.  Mech.  &  Engng., 
80,  59-64. 

Durlofsky,  L.J.  Engquist,  B.  &  Osher,  S.J.  1990  "Triangle-based 
TVD  schemes  for  hyperbolic  conservation  laws",  Unpublished  CAM 
Report  89-33. 

E.  W.  &  Hou,  T.Y.  mo  Comm.  Pure  Appl.  Math.,  43,  821-855. 
Engquist,  B.  1987  In  Proc.  Conf.  on  Hyperbolic  Partial  Differential 
Equations,  (Eds.  C.  Carasso,  P.-A.  Raviart  &  D.  Serre).  pplO-22. 
Springer. 

Engquist,  B.  1989  7.  Comp.  Math.  7,  113-120. 

Engquist,  B.,  Greenbaum,  A.  &  Murphy,  W.D.  1989a  IEEE  Trans 
Magn.  25,  2804-2806. 

Engquist,  B,  &  Gustafsson,  B,  1987  Math.  Comp.  49,  39-64. 

Engquist,  B.  &  Halpern,  L.  1988  7.  Appl.  Math.  4,  21-45. 

Engquist,  B.  &  Halpern,  L.  1990  Math.  Meth.  Appl.  Sci.  13,  189- 
203. 

Engquist,  B.  Harten,  A.  &  Osher,  S.  1987  In  Proc.  Conf.  on  Large 
Scale  Scientific  Computing,  (Ed.  M.  Birkhauser)  ppl 97-208. 

Engquist,  B.  &  Hou  T.Y.  1988  In  Vortex  Methods,  (Ed.  C.R. 
Anderson  &  C.  Greengard).  pp68-82.  Springer. 

Engquist,  B.  &  Hou  T.Y.1989  SIAM  J.  Hum.  Anal.  26,  289-319. 
Engquist,  B.  &  Huynh,  Q.Q.  1990  In  Proc.  Numer.  Analysis  Conf., 
Merida  (Ed.  S.  Gomez),  pp  60-74. 

Engquist,  B.,  Lotstedt,  P.  &  Sjogreen,  B.  1989b  Math,  of  Comp.  52, 
509-537. 

Engquist,  B.  &  Sjogreen,  B.  1990  "Numerical  approximation  for 
hyperbolic  conservation  laws  with  stiff  terms".  Unpublished  CAM 
Report  89-07. 

Fatemi,  E.,  Gardner,  C.L.,  Jerome,  J,,  Osher,  S.J.  &  Rose,  D.J.  1991a  In 
Proceedings  of  the  Urbana  Workshop  on  Computational  Electronics, 
(Eds.  K.  Hess,  J.P.  LeBurton,  U.  Raviolo).  Kluver. 

Fatemi,  E.,  Jerome,  J.,  &  Osher,  S.J.  1991b  IEEE,  10,  232-243. 

Haeri,  M.B.,  Putierman,  S.J.,  Garcia,  A.  &  Roberts,  P.H.  1992  Phys. 
Rev.  Lett,  submitted. 

Halpern,  L.  1991,  SIAM  J.  Math.  Anal,  22,  1256-1283. 

Harten  A.  1988  7.  Comp.  Phys.  83,  148-184. 

Harten  A.  1990  In  The  Numerical  Modelling  of  Nonlinear  Stellar 
Puh  itions,  (Ed.  R  Buchler),  NATO  ASI  Series  C,  302,  239-262. 


Harten  A.  &  Chakravarthy,  S.R.  1992  J.  Comp.  Phys.,  submitted. 
Harten,  A.,  Engquist,  B.,  Osher,  S.  &  Chakrovarthy,  S.  1987  J. 
Comp.  Phys.  71^  231-303. 

Hills,  R.N.,  Loper,  D.E.  &  Roberts,  P.H.  1987c  Quart.  J.  Mech.  &  Appl. 
Math.  40,  279-301. 

Hills,  R.N.  &  Roberts,  P.H.  1987a  J.  Nonequilib.  Thermodynam.  12, 
165-181. 

Hills,  R.N.  &  Roberts,  P.H.  1987b  J.  Nonequilib.  Thermodynam.  12, 
183-195. 

Hills,  R.N.  &  Roberts,  P.H.  1988a  Int.  Comm.  Heat  Mass  Transfer.  IS, 
113-119. 

Hills,  R.N.  &  Roberts,  P.H.  1988b  J.  Engng  Math.  22,  93-106. 

Hills,  R.N.  &  Roberts,  P.H. 1988c  Int.  J.  Nonlin.  Mech.,  23,  327-339. 

Hills,  R.N.  &  Roberts,  P.H.  1989  J.  Low  Temp.  Phys.  74,  185-203. 
Hills,  R.N.  &  Roberts,  P.H.  1990  Int.  J.  Nonlin.  Mech.,  25,  319-329. 

Hills,  R.N.  &  Roberts,  P.H.  1991  Stab.  &  Appl.  Anal.  Cont.  Media,  1, 

1-26. 

Hills,  R.N.  &  Roberts,  P.H.  1991  Stab.  &  Appl.  Anal.  Cont.  Media,  1 , 
205-212. 

Jones,  C.A.,  Putterman,  S.J.  &.  Roberts,  P.H.  1986  7.  Phys  A19,  2991- 
3011. 

Jones,  C.A.,  &  Roberts,  P.H.  1990  Geophys.  Asirophys.  Fluid  Dynam., 
55, 1246-268. 

Katzer,  E.  &  Osher,  S.  1988  "Efficient  implementation  of  essentially 
non-oscillatory  schemes  for  systems  of  nonlinear  hyperbolic 
differential  equation;.  Unpublished  CAM  Report  88-14. 

Kreiss,  H.-O.  1992  Acta  Numerica,  to  appear. 

Kreiss,  H.-O.  &  Abrahamson,  L.  1992  "Embedded  Runge-Kuita 
methods  for  oscillatory  systems",  BIT,  to  appear. 

Kreiss,  H.-O.  &  Lorenz,  J.  1992,  "Manifolds  of  slow  solutions  for 
highly  oscillatory  problems",  Indiana  U.  Math.  J.,  submitted. 

Kreiss,  H.-O.  ,  Lorenz,  J.  &  Naughton,  M.J.  1991,  Advan.  Appl.  Math., 
12,  187-214. 

Kreiss,  H.-O.  &  Scherer,  G.  1992  "Methods  of  lines  for  hyperbolic 
differential  equations",  SIAM  J.  Numer.  Anal.,  to  appear 
Kreiss,  H.-O.  &  Wu,  L.  1992  to  appear. 

Kuang,  W.  &  Roberts,  P.H.  1990  Geophys.  Asirophys.  Fluid  Dynam., 
55, 199-239. 

Kuang,  W.  &  Roberts,  P.H.  1991  Geophys.  Asirophys.  Fluid  Dynam., 
60,  295-333. 

Kuang,  W.  &  Roberts,  P.H.  1992  Geophys.  Asirophys.  Fluid  Dynam., 
accepted. 

Lafon,  F.  &  Osher,  S.J.  1991  J.  Comp.  Phys.  96,  110-142. 


Lange,  C.G.  &  Smith  D.R.  1988  Studies  in  Appl  Math.  79,  1-63. 

Lange,  C.G.  &  Smith  D.R.  1988  Studies  in  Appl  Math,  submitted. 

Lange,  C.G.  &  Miura,  R.M.  1991  Studies  in  Appl  Math.  84,  231-273. 
Lange,  C.G.  &  Weinitschke  1991a  Studies  in  Appl  Math.  82,  7-41. 

Lange,  C.G.  &  Weinitschke  1991b  Math.  Mech  71,  T751-T753. 

Larraza,  A.  1987  "Universal  power  spectra  for  wave  turbulence: 
applications  to  wind  waves,  flicker  noise,  solar  wind  spectrum,  and 
classical  second  sound”.  Unpublished  CAM  Report  87-14, 

Larraza,  A.  Garrett,  S.L.  &  Putterman,  S.  1990  Phvs  Rev.  A41,  3144- 
3155. 

Lin  R-Q.  1990  Geophys.  Astrophvs.  Fluid  Dynam.,  54,  161-188. 

Liu,  X.-D.  1991  "A  maximum  principle  satisfying  modification  of 
triangle-based  adaptive  stencils  for  the  solution  of  hyperbolic 
conservation  laws".  Unpublished  CAM  Report  91-11 
Loper,  D.E.  «fe  Roberts,  P.H.  1987a  In  Structure  &  Dynamics  of 
Partially  Solidified  Systems  (Ed.  D.E.  Loper),  pp. 293-323.  Reidel. 
Loper,  D.E.  &  Roberts,  P.H.  1987b  in  Hydro  dynamical  Behavior  <&. 
Interacting  Particle  Systems,  (Ed.  G.  Papanicolaou),  pp.  11 3-1  16, 
Springer. 

Mao,  De-Kang  1988a  "A  treatment  of  discontinuity  in  shock 

capturing  finite  difference",  Unpublished  CAM  Report  88-08. 

Mao,  De-Kang  1988b  "A  treatment  of  discontinuity  in  shock 

capturing  finite  differences  I.  Single  conservation  laws".  Unpublished 
CAM  Report  88-32. 

Mao,  De-Kang  1990  J.  Comp.  Phys.,  92,  422. 

Mao,  De-Kang  1991  "A  treatment  for  discontinuities  tor  finite 

difference  methods",  J.  Comp.  Phys.,  accepted. 

Mao,  De-Kang  1990  "A  treatment  for  discontinuities  for  finite 

differences  methods  in  two-dimensional  case".  Unpublished  CAM 
Report  90-30. 

Marquina  A.  1992  SIAM  J.  Scient  &.  Statist.  Comp.,  accepted. 
Marquina  A.  1993  In  Proceedings  of  4th  International  Conference 
on  Hyperbolic  Problems:  Theory,  Numerical  Methods  &  Applications, 
accepted. 

Mostrel  M.  M.  1988  In  Proceedings  of  ASME  Symposium  on 
Advances  in  Fluid  Dynamics  ,  (Ed.  O.  Baysal). 

Mostrel  M.  M.  \9S9  Math,  of  Comp.,  52,  587-613. 

Mulder,  W.  1989a  7.  Comp.  Phys.  83,  303-323. 

Mulder,  W.  1989b  In  Proceedings  of  the  7/'^  International 
Conference  on  Numerical  Methods  in  Fluid  Dynamics,  Lecture  Notes 
in  Physics,  323.  Springer:  Heidelberg,  429-433. 

Mulder,  W.  1989c  In  Proceedings  of  the  Fourth  Copper  Mountain 
Conference  on  Multigrid  Methods,  (Eds.  J.  Mandel,  S.F.  McCormick,  J.E. 


Dendy,  C.  Harhat,  G.  Lonsdale,  S.V.  Farter,  J.W.  Ruge  &  K.  Siuben.) 
SIAM  Publication,  pp  348-364. 

Mulder,  \V,  1989d  In  Proceedings  of  AIAA  Conference  on 

Computational  Fluid  Dynamics,  Paper  89-194  . 

Mulder,  W.  1992  "A  high-resolution  Euler  solver",  J.  Comp.  Phys.. 
in  press. 

Mulder,  W.,  Osher,  S.  &  Sethian,  J.A.  1992,  "Computing  interface 
motion  in  compressible  gas  dynamics",  J.  Comp.  Phys.,  submitted. 

Ong,  E.G.  1992  "Uniform  refinement  of  a  tetrahedron”,  SIAM  J. 
Scient.  &  Statist.  Comp.,  to  appear. 

Osher,  S.J.  1991  In  Computing  Optimal  Geometries,  (Ed.  J.E.  Taylor). 
pp35-39.  AMS  Publication. 

Osher,  S.J.  &  Shu,  C.W.  1988  In  Lecture  Notes  in  Engineering  43. 
(Eds.  C.  C.  Chao,  S.A.  Orszag  &  W.  Shyy)  pp31 6-333.  Springer. 

Osher,  S.J.  &  Shu,  C.W.  1991  J.  Numer.  Anal.  28,  907-922. 
Papageorgiou,  D.T.  &  Smyrlis,  Y.S.  1991  "The  route  to  chaos  for  the 
Kuramoto-Sivashinsky  Equation",  Unpublished  CAM  Report  90-26. 
Putterman,  S.  &  Roberts,  P.H.  1988  Phys.  Repts.  168,  209-263. 
Putterman,  S.,  Roberts,  P.H.  &  Fiszdon,  \V.  1988,  Ph\s.  Lett. 
A 128,  203-206. 

Putti,  M.,  Yeh,  W.-G.,  &  Mulder,  VV.,  1992,  "A  triangular  finite 
volume  approach  with  high  resolution  upwind  terms  for  the  solution 
of  groundwater  transport  equations",  J.  Comp.  Phys.,  submitted. 

Rizzi,  A.  &  Engquist,  B.  1987  Math.  Comp.  72,  1-69. 

Roberts,  P.H.  1987  "Nonlinear  Dynamics  of  Rotating  Magnetic 
Systems".  Unpublished  CAM  Report  87-23. 

Roberts,  P.H.  1991  in  The  Sun  Cool  Stars;  Activity,  Magnetism, 
Dynamos,  (ed.  I.  Tuominen,  D.  Moss  &  G.  Rudiger)  pp37-56.  Springer. 
Roberts,  P.H.  &  Loper  D.E.  1987  in  Structure  <Sc  Dynamics  of 
Partially  Solidified  Systems  (Ed.  D.E.  Loper),  pp229-290.  Reidel. 

Shu,  C.-W.  Erlebacher,  G.,  Zang,  T.A.,  Whitaker,  T.,  &  Osher,  S.  1992, 
submitted. 

Shu,  C.-W.  8l  Osher,  S.J.  1989  J.  Comp.  Phys.  83,  32-78. 

Sjogreen,  B.  1988  "Shock  Computations  on  Parallel  Computers", 
Unpublished  CAM  Report  87-27. 

Ton,  V.T.,  Karagozian,  A.R.,  Engquist,  B.E.  &  Osher,  S.J.  1992 
"Numerical  simulation  of  inviscid  detonation  waves  with  finite  rate 
chemistry".  Unpublished  CAM  Report  91-22. 

Wu,  L.  1992  "Stability  of  difference  approximations  for  initial 
boundary  value  problems",  submitted. 

Yang,  H.  1990  J.  Comp.  Phys.,  89,  125-160. 


